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Abstract 

An analysis of the CP violating asymmetry in — > (3-7r) ± decays in 
the Standard Model and, by means of the mass insertion approximation, 
in a wide class of possible super symmetric extensions, is presented. We 
find that the natural order of magnitude for this asymmetry is 0(1O -5 ) in 
both cases. Within super symmetric models effects as large as 0(1O -4 ) are 
possible, but only in a restricted range of the relevant parameters. 
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1. The origin of CP violation is one of the fundamental questions of particle 
physics and cosmology which remains an open problem to date. The recent 
measurements of e'/e |]J (see also 0) represent an important step forward in our 
understanding of this phenomenon, since they have ruled out superweak scenarios. 
Nonetheless we are still far from a quantitative description of the dynamics which 
generate the amount of CP violation observed in hadronic processes. Indeed, even 
within the Standard Model (SM) it is very hard to predict the value of e'/e in 
terms of the parameters of the Cabibbo-Kobayashi-Maskawa (CKM) matrix (see 
e.g. Refs. and references therein). Given the large theoretical uncertainties 
affecting the calculation of this quantity, it is very useful to collect additional 
experimental information about CP violation in \AS\ = 1 transitions. In this 
respect charge asymmetries in non-leptonic decays, as the difference in K + — > 
(37r) + and K~ — > (37r) _ Dalitz plot distributions ||, represent an interesting 
class of observables, since they are straight direct CP-violating effects free from 
| AS | = 2 contaminations. Moreover, contrary to e'/e, these asymmetries stem 
from the interference of two AI = 1/2 amplitudes and do not necessarily suffer 
the suppression of AI = 3/2 transitions. In spite of these advantages, however, 
within the SM such observables are expected to be very small, of (9(10 -5 ), due 
to the constraints from e'/e and the smallness of final-state interactions |]. IJ. 
A natural question is whether extensions of the SM could enhance these CP- 
violating asymmetries at such a level that they could be recognized as a clear 
signal of new physics. 

Good candidates to provide new large CP violating effects are the supersym- 
metric extensions of the SM with generic flavour couplings and minimal particle 
content. In this framework, among the possible contributions which may be en- 
visaged, it has been recently recognized the importance of the chromomagnetic 
operator (CMO). Its CP-odd contribution can become large in the presence of 
misalignment between quark and squark mass matrices, and, without conflict 
with the experimental determination of the K°-K° mixing amplitude, it can ac- 
count for the largest part of the measured e'/e f7|]-|PI|]. Actually a non-standard 
CMO is the only possibility, within this framework, to considerably affect the 
CP- violating part of AI = 1/2 amplitudes without serious fine-tuning problems 



in | AS 1 1 = 2 processes [10 



In this paper we investigate the possibility of using the CMO to enhance 
CP violating effects in K — > 3tt decays. We work under the assumption that 
the Wilson coefficient of this operator is mostly due to left-right mixing among 
down-like squarks. This implies that a large chromomagnetic term is necessarily 
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accompanied by sizeable corrections to e and Ki — > 7r°e + e~. We thus perform 
a combined analysis of these processes together with e'/e and K — > 3tt decays. 
We have not considered other possible supersymmetric sources of CP violation 
in K — > 37T decays, such as left-left or right-right squark mixing, since only the 
left-right ones trigger the enhancement of the CMO which is the most promising 
candidate to give an observable effect. 

Our main conclusions are that, even within supersymmetry, effects at the 
level of 0(1O -4 ) may be observed only under special circumstances, among which 
large cancellations of different contributions in e'/e. Otherwise, similar to the 
SM, the natural order of magnitude of the charge asymmetries will remain of 
0(1O~ 5 ). As we shall discuss, this conclusion is rather general and applies also 
to other direct CP-violating observables in non-leptonic processes. 

The paper is organized as follows. We first derive the main formulae needed 
to evaluate the CMO contribution to CP-odd asymmetries in K ± — > {Sir^ de- 
cays. Then we discuss the role of down-type left-right mass insertions, including 
contributions from the CMO, in the K°-K° mixing amplitude. Finally a com- 
bined analysis of CP- violating effects in K ± — > (37r) ± decays, taking into account 
the constraints imposed by e, e'/e and Kl — > n e + e~, is presented. The results 
are summarized in the conclusions. 



2. We start by analyzing the charge asymmetries in —> 7r ± 7r ± 7r =F decays. 
As discussed in Ref. || ||, the most interesting CP- violating observable is the 
asymmetry in the Dalitz plot slopes g±. Neglecting the suppressed A J = 3/2 
contributions, this can be written as 



9+ -9- 
9+ + 9- 



\mb Ima 



Re6 Rea 



sin(a - A)) 



(1) 



where the weak amplitudes a and b are defined by the momentum expansion of 
A(K + — > 7r + 7r + 7r~) around the center of the Dalitz Plot, 



A(K~' 
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and ao, A) are the small rescattering phases, known from chiral perturbation 
theory (ChPT) [ill], evaluated at Y — 0. In the limit where we neglect AI = 3/2 



2 



contributions, the slope asymmetries of K — > 7r ± 7r ± 7r T and K — > ti ^ ^ 
modes are identical. 

Since a and b are AI = 1/2 amplitudes, on general grounds one expects 
Ima/Rea and Imb/Keb to be both of the same order as the weak phase of Aq = 
A(K — > (27r)/=o), namely 

Ima Im& Imv4 ^ ( e'\ t _ t R,eA in _4 



Re — ~ 10~ 4 . (4) 

Rea Refe ReA U / ReA 2 



Given that sin(«o — Po) ~ 0.1 PHI , this sets the "natural" order of magnitude 
for the asymmetry to 10~ 5 ||. Actually, within the SM, the situation is even 
worse: neglecting the CMO, which in this case has a very small coefficient, the 
asymmetry vanishes at the lowest order in the chiral expansion. This happens 
because at this order there is only one octet operator which generate the same 
weak phase to all the AI = 1/2 amplitudes . Clearly the situation may improve 
if the contribution of the CMO is enhanced by supersymmetric effects, which we 
now discuss. 

Let us start from the relevant piece of the effective Hamiltonian. This can be 

Hmag = C+Q+ + C-Q- + h.C. , (5) 



written as fTO 
where 



Qt = (s L a^t a G%d R ± s R a^t a G%d L ) (6) 



9_ 

167T 2 

and the dominant contribution to Wilson coefficients, generated by gluino ex- 
change diagrams, is given by [|K| [L2| 



± = na s{ m~ 9 ) 



K«) 21 ± G (x gq ) . (7) 



Here (S^ R )ij = (My) iL j R /m 2 ~ denote the off-diagonal entries of the (down-type) 
squark mass matrix in the super-CKM basis |13j and x gq = m|/m~ the ratio of 
gluino and (average) squark mass squared. The explicit expression of Gq(x) can 
be found in [Hj. 

The realization of in terms of meson fields, to the lowest order in 1/N C 
and in the derivative expansion, can be written as 



, 11 PM 2 k 

of- 



/g 2567T 2 m s + m<i 
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where U = exp(i20// 7r ), f n = 132 MeV, <fi is the octet field of pseudoscalar 
mesons, and the overall coupling has been fixed by the chiral quark model estimate 
of Ref. 



I4| . Using Eq. ([8]) we can derive the following matrix elements 
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(9) 
(10) 

(11) 



The .B-factors, i? si , have been introduced to parametrize our ignorance of the pre- 
cise overall coefficient in Eq. (||) and of possible higher-order terms. For practical 
purposes, in the following we shall also set m s + m d = 110 MeV in Eqs. (P|)-(|TT|), 
encoding in the B g i the remaining uncertainty on the true value of the quark 
masses. [] 

For the K — > 3tc amplitudes defined in (0) we obtain 



Ima 



Rea 
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2M|G 8 16vr 2 /2(m s + m d ) 



B g3 lmC, 



+ 



Im6 

Reb 



0, (12) 



where we have used the lowest-order chiral relation between Rea and Rev4 , ex- 
pressing the latter in terms of the standard coupling G 8 = 9.1 x 10~ 6 GeV -2 (i.e. 
Rea = 2Mj£Grg/3). In view of the numerical analysis, it is convenient to introduce 
the following simple expression, which can be readily derived from Eq. ([T2] ) and 



^mag in 



(1) 



9+ - 9- 


~ 1.97 x 


9+ + 9- 





rj a s (rrig) 500GeV G (x 



<!<! ) 



a s (500GeV) m § 

= (6E R ) 21 ± (8E R y l2 



\B g3 lm5^ 



(13) 



where we have defined 8± = {5^ R ) 2 i ± (^L/?)i2 = (^ R ) 2 i ± {5 RL ) 2 \- We found 
very useful to introduce 5± since these are the natural couplings appearing at 
first order in any parity conserving (+) or parity violating (— ) observable. In the 
evaluation of the numerical coefficient above we have used a s (500GeV) = 0.096 
and, as anticipated, m s + m d = 110 MeV. The parameter rj ~ 0.9 [TIJ is the 
correcting factor due to the running of the Wilson coefficient from m g to the 
operator renormalization scale. 

1 For an extensive discussion about the possible chiral realizations of the CMO see Ref. Jll| . 
The factor B g z could in principle be a function of Y, but for simplicity in the following we will 
ignore this possibility; B g i coincides with the Bq of (l(J for m s + nid = 110 MeV. 
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3. An important constraint on the couplings 5± comes from K°-K° mixing. 
Besides the usual gluino-box amplitudes widely discussed in the literature (see 



e.g. Refs. [12|, 16, |H| and references therein) further contributions arise from the 
single and double insertion of the CMO. Schematically we can write 



A susy (K° -> K c 



A 



boxes 



+ -4.1-mag + A2- 



2— mag • 



(14) 



vAboxes, which is dominated by short-distance contributions, can conveniently be 
written in the form 



A 



boxes 



M 



K 



m| 432 \m s + rrid 



M 2 K f 2 K xU 6 (x gq ) (85 V2 B 2 + 3 V3 B 3 ) (5 2 + + Si) 



+ x gq h{x gq ) (33^54 + 15^5^5) [8\ - 5 2 _) 



(15) 



where the definitions of fe,(x), f§{x) and of the Bi parameters have been taken 

lalysis we will use 

a s (m~ g ) 500GeV" 



from Ref. [IB|. In the numerical analysis we will use the simplified expression 
Aoxes = 2.9 x KT n GeV 2 x 



a,(500GeV) m- g 



B + 51 + B_ 5t 



(16) 



where B± are coefficients of 0(1), which may vary by a factor of 2 3 depending 
on the values of the B parameters of the AS" = 2 operators, the precise value of 



x gg , the perturbative QCD corrections, etc. [16 



The single and double insertions of the CMO are expected to be dominated 
by long distance contributions. For illustration we give here the expression of the 



single 7r u contribution, already considered in Ref. [19 



At 
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Ml 
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K 



Ml 



Ml I 32vr m s + m d 



G 8 f 2 ^ v G (x gq )5 + . (17) 



m g 



One can then generalize the above expression to include the contribution from 
other one-meson states, such as the r\ and rj'. In the following we will use the 
simplified expression 



Ai. 



mag 



4.8 x KT 13 GeV 2 x 



i] a s (rrig) 500GeV G (x gq ) 
a s (500GeV) m g G (l) 



;is) 
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where the numerical coefficient has been computed from Eq. (|T7D and we have 
absorbed the hadronic uncertainties, namely B g i and contributions from inter- 
mediate states other than the ir°, in the factor K\. Using the Gell-Mann-Okubo 
mass formula the tt° and rj contributions would cancel, thus one may argue that 
K\ should be substantially smaller than one. We know, however, that a similar 
argument fails for Kl —>■ 77, where the effective coupling, corresponding to our 
Ki, is of 0(1). Note that intermediate states with parity opposite to the one 
pseudoscalar-meson state may give contributions proportional to <5_. We have 
neglected these effects in our analysis. 

Similarly, in the case of the double insertion one gets 



2— mag 



(^ |ftmagk°> 



Ml - Ml 



(7T°|H mag |^ ) 



(19) 



that proceeding as before leads to 



A 



2— mag 



1.9 x KT^GeV 2 x 



n 2 
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k 2 S 2 + . (20) 



4. For the other two quantities which are used in our analysis, namely e'/e and 
the BK(Ki 7r°e + e - ), rather than giving the explicit analytic expressions, for 
which we refer the reader to Ref. [10|, we only list here two convenient expressions: 

r] a s (rrig) 500GeV Go(x gq ) 



Re 



e 



92.6 x 



mag 



BR(K L 



7r°e + e- 



/ mag 
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^ 7 (500GeV, l)Go(l) 



B g2 Im5_ 



B|(Im5 H 



(21) 



(22) 



7r°e + e 



where the definitions of Bt and can be found in Ref. PH[ , and in Kl 
we have neglected the interference with the SM contribution. 

Besides the numerical expressions given above, for our study we also use the 
following experimental inputs: 



(2.28 ±0.02) x 10 



-3 



Re It 



(21.2 ±4.6) x 10 



-1 



BR(K L 



7T°e + e 



< 5.6 x 10 



-10 



[18] 



(23) 
(24) 
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We are now ready to discuss the bounds in the Im5 + -Im5_ plane imposed by 
the experimental measurements and the theoretical expressions given in Eqs. (|16D . 
(|IBD , (p0|)-(|2"2|). As usual, we impose the constraints by requiring that all CP vio- 
lating observables are saturated by the supersymmetric effects considered above, 
i.e. neglecting the SM contributions. By setting the values of all hadronic pa- 
rameters (namely the B^s and the k^s) to one, at our reference values of gluino 
and squark masses {m g = 500 GeV and x gg = 1), we find Im<5_ ~ lm5 + ~ 10~ 5 . 
According to Eq. (0), this implies that the K — > 3tt asymmetry is of the same 
order, as found within the SM. It is interesting to note that values of Im^ ~ 10~ 5 
are consistent with the approximate flavor-symmetry scenario of Ref. J7J, where 
\5±\ ^ sm6 c m s /mq. In this framework one could therefore find a "natural" su- 
persymmetric explanation for both e and e' . 

To obtain larger values of the charge asymmetry in K — ► 3tt decays one 
has to relax the bound on Im<5 + . To this purpose we note that e'/e put an 
explicit constraint only on Im5_ but not on Im<5 + , whereas BR,(Kl — > 7r°e + e~) 
and ^4boxes put upper bounds on Im<5 + only at the level of 10~ 4 . The strongest 
limit is set by the contribution of *4.i_ ma g to e, since this term is linear in lm5 + (the 
quadratic terms may become competitive only for KeS/lm8 ^> 1 or if K\ <C k?). 
This contribution is subject to a large uncertainty which is parametrized by K\. 
Therefore one can relax the bound by taking for K\ a value sensibly smaller than 
one, as for example done in Ref. [JHJ. In Fig. ffl, we display the bounds obtained 
for Ki = 1.0, 0.3 and 0.1. Only in the latter (optimistic) case one may obtain 
K — > 3tt asymmetries in the 10~ 4 range. This is very similar to what has been 
found in Ref. |19[ for the CP asymmetries of hyperon decays. 

Even accepting values of Im<5 + ~ 10~ 4 , it remains to be explained the large 
cancellation between lm(S^ R )2i and Im(5^ L ) 2 i necessary to satisfy the 10~ 5 bound 
on ImS- imposed by e'/e. An underlying mechanism forcing this cancellation 
exists, however, in the U(2) models considered in Ref. 0. 

In principle one could relax the e'/e constraint on Im£_ by allowing other 
contributions to e'/e to be large. For example if one accept the striking result 
of Ref. |2(J,0 (£'/s)sm ~ —120 x 10~ 4 , one would need a large supersymmetric 



contribution, corresponding to Im5_ ~ 10 4 to reconcile the theory with the 
experimental number. This, however, leads to a new fine tuning problem, because 



2 This value completely disagrees in sign and size with experimental measurements and with 
many theoretical estimates. We consider indeed premature to use it in phenomenological anal- 
yses. 
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Figure 1: Constraints in the Im5 + -Im5_ plane imposed by e'/e, e and 
BR(Kl — > 7T°e + e~). All bounds are obtained for x gq = 1 and scale linearly 
with (500GeV/mg). The vertical dotted lines correspond to the e constraint on 
^i-mag for Ki = 1.0, 0.3 and 0.1 (from left to right); the other vertical lines are 
obtained from the e bound on (.Aboxes + -42-mag) for 5- = 0, Re<5 + = Im<5 + and 
B+ — K2 — 1 (dashed) or B + — — k% — 1 (dot-dashed). The constraint from 
BR,(Kl — > 7r°e + e _ ) is obtained for Bt = 1. The limits from e'/e are obtained 
for B g2 = 1 and (e'/e)sM — (horizontal dash-dotted line) or the extreme case 
(£'/ £ )sm < —120 x 10~ 4 pOf (horizontal shadowed region). 
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then the natural order of magnitude of e'/e would be 10 2 rather than 10 3 . 

5. Large CP violating effects triggered by misalignments between quark and 
squark mass matrices are among the most promising phenomena to uncover Su- 
persymmetry at low energy. Among the possible effects those driven by the 
chromomagnetic operator are particularly interesting since they could completely 
account for the already measured CP violating parameters, e and e'/e. In this 
paper we have studied the possibility that effects of the CMO are detectable from 
the enhanced asymmetry in K — > 3tt decays. We find that this is possible only 
if several conditions, on which we have a poor theoretical control, conspire in 
the same direction. Moreover, even if this were the case, one should then face a 
fine-tuning problem in e'/e. 

Our analysis of supersymmetric CP-violating effects in K — > 37r decays is 
parallel to those recently performed in K — > tttt^ |2l] and hyperon decays ||19|| . In 



all these cases the conclusions are hampered by poor knowledge of some hadronic 
parameters, which in the future will hopefully be computed on the lattice. We 
stress that this problem is absent (or at least much simpler) in rare K decays like 
Ki — > 7r°z/z/(e + e - ) [IIJ, whose experimental investigation will definitely provide 



useful and unambiguous information about the nature of CP violation. 

We thank I. Mannelli for interesting discussions that stimulated us to start 
this work. G.I. thanks LAL and the Laboratoire de Physique Theorique of the 
CNRS at Universite de Paris XI for hospitality during the completion of this 
work. 
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